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MODAL  THEORY  OF  TRANSVERSE  ACOUSTIC  COHERENCE  IN  SHALLOW 

OCEANS 


1.  INTRODUCTION 

Ocean-acoustic  systems  operating  in  continental-shelf  regions  must  cope  with  dynamic  sound-speed 
fields  that  are  driven  by  oceanographic  processes.  In  such  environments,  acoustic  fields  are  known  to  fluctu¬ 
ate  significantly,  and  it  is  clear  from  the  physics  of  the  problem  that  the  acoustical  fluctuations  are  governed 
by  those  of  the  underlying  sound-speed  field,  which  is  affected  by  diverse  oceanic  processes  including  sur¬ 
face  waves,  internal  waves,  bubble  clouds,  fish  schools,  bottom  topography,  and  sub-bottom  geology.  The 
present  work  is  focused  solely  on  the  relationship  between  the  internal-wave  spectrum  and  acoustic  coher¬ 
ence  in  shallow  oceans. 

Awareness  of  the  impact  of  internal  waves  and  tides  on  acoustic  propagation  in  the  ocean  was  attributed 
by  Whitmarsh  et  al.  [1]  to  work  conducted  during  World  War  II.  Whitmarsh  et  al.  themselves  compared 
acoustic  scattering  of  25  kHz  pulses  at  256  m  range  with  theoretical  predictions  of  Skudrzyk  [2]  that  were 
based  on  spatial  correlation  functions  and  the  focusing  of  sound  by  anomalous  sound-speed  patches.  Their 
conclusions,  valid  for  high-frequency  propagation  at  short  range,  are  that  amplitude  variance  increases  with 
range,  independent  of  the  sea  state.  Lee  [3,  4]  based  his  efforts  on  work  done  by  Ufford  that  was  also 
done  during  World  War  II  and  published  afterward  [5,  6].  Using  ray-based  numerical  simulations,  Lee  [3] 
demonstrated  that  internal  waves  can  cause  an  increase  in  acoustic  intensity  contrast  in  comparison  with  the 
absence  of  such  waves.  DeFerrari  [7]  conducted  ray-based  simulations  of  long-range  deep-water  acoustic 
propagation  through  internal  tides  and  waves  that  demonstrated  phase  shifts  among  multipath  arrivals  are 
the  underlying  cause  of  temporal  fluctuations  in  acoustic  intensity. 

Breaking  with  the  ray-based  work,  much  of  it  in  the  temporal  domain,  Flatte  and  Tappert  [8]  combined  a 
parabolic -equation  acoustic  model  with  the  Garrett-Munk  internal-wave  spectrum,  in  a  study  which  demon¬ 
strated  significant  impact  of  internal  waves  on  acoustic  propagation  in  deep-water  environments.  Munk  and 
Zachariasen  [9]  combined  a  geometric-optics  approximation  of  the  acoustic  field  with  the  Garrett-Munk 
spectrum  to  rigorously  relate,  for  the  first  time,  internal-wave  and  acoustical  statistics,  which  compared  well 
with  data  from  two  experiments.  A  good  historical  review  of  work  in  this  area  may  be  found  in  Flatte  [10]. 

In  ocean  basins  with  a  deep  sound  channel,  the  physical  difficulties  related  to  boundary  interactions 
are  somewhat  less  daunting  than  in  shallow-water  environments,  which  have  multiple  boundary  reflections. 
Here,  path-integral  techniques  were  successfully  adapted  to  the  propagation  of  waves  in  random  media  by 
Dashen  [11],  Using  this  approach,  Esswein  and  Flatte  rigorously  related  sound-speed  fluctuations  due  to 
linear  internal  waves,  to  the  spatial  and  temporal  coherence  of  acoustic  fields  [12]. 

Path-integral  theories  continue  to  be  refined  and  used  today  for  deep-water  work  being  done  by  the 
Acoustic  Thermometry  of  Ocean  Climate  (ATOC)  and  North  Pacific  Acoustic  Laboratory  (NPAL)  groups 
[13-15].  One  significant  finding  of  the  ATOC  effort  [16]  has  been  the  experimental  verification  in  deep 
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water  of  the  path-integral  prediction  that  the  mutual-coherence  function  of  the  acoustic  field  with  complex 
pressure  P(yi)  has  the  form 

(P1P2)  =  (|JPi||P2|)exp(-iD(jI -y2)),  (1) 


where  the  phase-structure  function  D{y\  —  y2)  is  the  variance  of  stochastic  phase  differences  between  two 
points  in  the  acoustic  field,  and  (\P\  1 1 /-*2 1 )  is  the  amplitude  covariance.  One  objective  of  the  present  work  is 
to  determine  whether  this  exponential  form  holds  in  shallow  water. 

A  normal-mode  approach  to  acoustical  fluctuations  has  also  received  long-standing  attention,  with  the 
general  conclusion  that  mode  coupling  induced  by  internal-wave  fluctuations  plays  an  important  role  in  both 
deep  and  shallow  water.  Such  couplings  were  clearly  illustrated  by  the  parabolic  equation  simulations  of 
Colosi  and  Flatte  [17].  Although  it  is  an  important  component  of  the  problem  that  eventually  needs  atten¬ 
tion,  the  inclusion  of  coupling  leads  to  considerable  complications  that  require  numerical  evaluation  [18]. 
By  neglecting  mode  coupling,  the  analysis  can  be  completed  analytically,  thus  providing  physical  insight 
into  the  relationship  between  acoustical  and  environmental  fluctuations.  In  this  spirit,  the  incorporation  of 
coupled  modes  into  the  path-integral  approach  is  not  included  in  the  present  work. 

With  respect  to  shallow-water  oceans,  various  other  parts  of  the  problem  have  received  attention  over 
the  past  fifteen  years.  Preisig  and  Duda  illustrated  the  resonant  interaction  of  acoustic  waves  as  they  pass 
through  large,  nonlinear  internal  waves  [19-21].  Finette  et  al.  [22]  performed  data-constrained  numerical 
simulations  of  acoustic  intensity  fluctuations  observed  in  the  1995  Shallow-Water  Acoustics  in  a  Random 
Medium  (SWARM)  experiment  [23].  Using  numerical  simulations,  Finette  and  Oba  [24]  discovered  the 
phenomenon  of  horizontal  refraction  of  acoustic  waves  propagating  nearly  parallel  to  nonlinear  internal- 
wave  crests.  However,  a  quantitative  prediction  of  the  transverse  mutual-coherence  function  remains  elusive 
for  shallow  oceans.  The  present  work  represents  a  step  in  that  direction. 

As  currently  formulated,  path-integral  theory  applies  to  the  fluctuation  of  a  single  refracted  Fermat 
path  that  does  not  interact  with  boundaries.  Clearly  this  restriction  poses  difficulties  for  shallow-water 
acoustic  propagation,  which  not  only  has  multiple  boundary  interactions,  but  is  in  fact  modal  for  the  lowest 
frequencies.  To  circumvent  this  problem  in  the  present  work,  the  acoustic  propagation  problem  is  formulated 
using  adiabatic  normal  modes,  and  path-integral  theory  is  applied  to  the  propagation  of  mode  amplitudes  in 
the  horizontal  plane  where  there  are  no  boundary  interactions  —  all  paths  are  refracted. 

In  the  present  work,  the  acoustic  field  is  separated  into  adiabatic  normal  modes  where  a  parabolic  ap¬ 
proximation  of  the  Helmholtz  equation  governs  the  propagation  of  mode  amplitudes.  Here  the  horizontal 
wavenumbers  provide  stochastic  indices  of  refraction.  The  wavefunction  satisfying  the  parabolic  equation 
is  then  written  as  an  integral  over  all  paths  between  a  source  and  receiver,  and  the  second  moment  is  ob¬ 
tained  by  multiplying  and  ensemble  averaging  the  wavefunctions  of  two  receivers  separated  in  space.  After 
application  of  the  Markov  approximation  and  the  Wiener-Khintchine  theorem,  the  resulting  phase-structure 
function  emerges  as  an  integral  over  the  spectrum  of  transverse  horizontal  wavenumbers.  Perturbations  in 
the  wavenumbers  are  related  to  those  of  the  sound  speed  using  an  eigenvalue  perturbation  method  adopted 
from  quantum  mechanics.  The  potential  sound-speed  gradient  is  then  used  to  express  the  perturbations  as  a 
linear  combination  of  internal-wave  modes.  The  final  result  expresses  the  desired  phase-structure  function 
as  an  integral  over  the  transverse  internal-wave  spectrum,  which  must  be  either  measured  or  modeled.  A 
shallow-water  version  of  the  Garrett-Munk  spectrum  is  then  introduced  to  provide  one  possible  model. 
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The  remainder  of  this  work  begins  in  Ch.  2  with  theoretical  development  that  expresses  the  acoustic 
coherence  function  as  a  weighted  sum  over  mode-amplitude  coherence  functions.  In  Ch.  3,  path-integral 
theory  is  applied  to  fluctuations  of  the  mode-amplitude  coefficients  to  obtain  an  expression  for  the  spectral 
representation  of  the  phase-structure  function.  In  Ch.  4,  fluctuations  in  the  acoustic  horizontal  wavenumbers 
are  related  to  those  of  a  linear  internal-wave  field.  In  Ch.  5,  everything  is  pulled  together  to  produce  a  set  of 
numerical  simulations  that  illustrate  the  properties  of  transverse  coherence  under  this  theory. 


2.  MODAL  REPRESENTATION  OF  COHERENCE 
2.1  Mutual-Coherence  Function 

In  the  works  of  Esswein  and  Flatte  [12],  Flatte  et  al.  [25],  and  Flatte  and  Stoughton  [26],  the  scattering 
strength  is  defined  as 


d>  = 


rSL 


k(x(s))/J.(x(s))ds 


(2) 


and  the  phase-structure  function  is  defined  as 


D(xi,x2)  =  /  (j^  k(xi(s))n(xi(s))ds  -  k(x2(s))/i(x2(s))dj 


(3) 


where  k(x)  is  the  unperturbed  wavenumber  of  the  background  environment,  jl(x)  are  stochastic  perturbations 
to  the  index  of  refraction,  .7(5)  are  spatial  coordinates  of  a  ray  path,  s  is  arc  length  along  the  unperturbed  ray, 
and  Sl  is  the  total  arc  length.  Properties  of  the  phase-structure  function  are  discussed  in  these  works. 


In  the  limit  of  geometric  optics  [25],  the  total  phase-fluctuation  8(j)  accumulates  along  the  propagation 
path. 


rSL 

8(f)  =  k(x(s))/j.(x(s))ds. 

Jo 


(4) 


Under  these  conditions  the  scattering  strength  becomes  the  variance  of  phase  fluctuations  <t>  =  (502),  and 
likewise,  the  phase-structure  function  becomes  the  variance  of  fluctuations  in  phase  difference  [13,  26], 
D(x \,x2)  =  ((8<p\  —  S^2).  The  structure  function  is  nonnegative,  and  for  large  hydrophone  separations  it 

approaches  twice  the  scattering  strength,  ^(50 1  —  f  — >  24>. 

Under  further  assumptions  of  Gaussian  statistics  for  the  phase  fluctuations,  it  has  been  shown  [25]  that 
the  normalized  covariance  between  two  points  of  a  complex  stochastic  field  P(x)  may  be  expressed  as  an 
exponential  of  the  phase-structure  function, 


PoMPo  (*2) 


exp  [~\D{x i,x2)]  , 


(5) 


where  Pq  is  the  field  in  the  absence  of  fluctuations.  This  well-established  equation  has  been  derived  using 
several  different  methods  [10,  27].  Equation  (5)  can  be  derived  under  less  restrictive  conditions  using  the 
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path-integral  approach  [10,  11,  27].  In  particular-,  Dashen  [11]  has  shown  that  the  only  requirements  for 
the  validity  of  Eq.  (5)  are  that  the  medium  be  isotropic  and  the  narrow-angle  parabolic  approximation  hold. 
Thus  Gaussian  statistics  for  the  phase  fluctuations  are  not  strictly  necessary.  Colosi  et  al.  [16]  has  shown 
that  Eq.  (5)  may  be  expressed  in  terms  of  amplitude  and  phase  correlations  as 

(P1P2)  =  <N  |P2|)exp  [-\D(x ux2)\  •  (6) 


Under  conditions  of  statistical  homogeneity  [15],  the  mutual-coherence  function  Rp(Ay )  of  the  acoustic 
pressure  P(y)  is  defined  by 


RP(Ay)  = 


(P(y)P*(y  +  Ay))tJ 

(ipwi2)w 


(V) 


where  Ay  is  hydrophone  separation  transverse  to  the  acoustic  propagation  path,  the  ensemble  average  is 
both  over  blocks  of  time  assuming  the  ergodic  hypothesis,  and  over  the  sub-diagonal  elements  of  a  Toeplitz 
cross-spectral  density  matrix.  By  combining  Eqs.  (6)  and  (7),  the  relationship  between  the  mutual-coherence 
function  and  the  phase-structure  function  is 


Rp(Ay) 


(|P(y)||P(y  +  Ay)|)^ 


exp[-yP>(Ay)]  . 


(8) 


To  remove  the  amplitude  correlation  from  measured  acoustic  data  for  comparison  with  the  theory  developed 
below,  the  complex  hydrophone  pressures  P,  are  normalized  as  p  =  P//|p|.  Thus  the  normalized  mutual- 
coherence  function  is  related  to  the  phase-structure  function  by 

Pp(Ay)  =  (P(y)P*(y  +  Ay))ty 

=  exp  [— ^D(Ay)]  .  (9) 

Note,  with  this  definition  the  reference  coherence  is  e  12  =  0.6,  where  the  coherence  scale  length  L  is 
defined  to  be  Ay  =  L  such  that  D  =  1. 

By  making  a  small-aperture  approximation  [12,  26,  27],  the  form  of  the  phase- structure  function  was 
determined  to  obey  a  power  law 

DP(Ay)=^  ,  (10) 

where  L  is  the  transverse  coherence  length.  The  correct  value  to  use  for  n  is  not  clearly  resolved  in  the 
literature.  Cox  [28]  simply  assumed  n  =  1  in  his  work  that  predates  the  path-integral  approach,  and  found 
reasonable  agreement  with  his  data.  Carey  [29]  reviewed  various  works  that  assumed  n  =  1.0, 1.5, 2.0,  and 
then  presented  an  array-signal-gain  approach  that  is  solvable  for  those  three  values.  For  observations  of 
transverse  coherence  in  the  Strait  of  Korea,  Carey  et  al.  [30]  used  n  =  2.0.  Esswein  and  Flatte  [12]  evalu¬ 
ated  Dp(Ay)  numerically  but  did  not  report  an  exponent.  Dashen  et  al.  [27]  fit  a  phenomenological  phase- 
correlation  function  with  exponent  n  =  1.5  to  Esswein’s  numerical  results.  This  value  was  then  adopted  in 
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the  follow-on  works  of  Flatte  and  Stoughton  [26]  and  Flatte  and  Rovner  [14].  For  deep-water  data  obtained 
in  the  Pacific,  Andrew  et  al.  [15]  found  good  agreement  between  the  observed  mutual-coherence  function, 
detailed  numerical  calculations  using  the  technique  of  Flatte  and  Rovner  [14],  and  a  small-aperture  empiri¬ 
cal  exponential  with  n  =  1.5.  Yet  there  are  theoretical  indications  in  Dashen  [11]  and  in  the  results  presented 
below  that  n  =  2.0  for  small  separations  of  the  path-integral  solution.  Whatever  the  correct  value  of  n  may 
be  for  small  Ay,  the  exponential  approximation  must  break  down  at  large  separations  [26]  where  the  power 
law  vanishes,  whereas  the  coherence  decays  to  a  constant  value,  Rp(Ay)  — >  exp(— <3>). 

2.2  Modal  Expansion  of  Acoustic  Field 

The  problem  is  formulated  in  cylindrical  coordinates  u  =  (r,  6 ,  z)  beginning  with  the  frequency-domain 
Flelmholtz  equation 


1  d 

r  dr 


d2P  _..d 
+  r2d02  +  P^dz 


1  dP 

p(z)  dz 


+  k2(u)P 


-Qs 


S(r)8(z-zs) 
2k  r 


(11) 


which  describes  the  acoustic  pressure  P(u)  of  a  narrowband  source  with  angular  frequency  co  =  2nf  and 
wavenumber  k(u)  =  (o/c{ u),  where  c(u)  is  the  sound  speed,  and  p(z)  is  the  mean  density  stratification.  Here 
the  source  term  takes  the  form  5(u  —  uv)  =  8(r')8(z  —  zs)/2nr,  with  Qs  as  the  source  strength.  A  separation 
of  variables  [31]  then  leads  to  an  expansion  of  the  acoustic  pressure  in  terms  of  local  modes,  which  vary 
over  the  horizontal  plane  [32-34], 


M 


P(r,d,z )  =  Y,Am(r,d)xVm(r,d;z). 

m—  1 


(12) 


The  real  eigenfunctions  'F „,{r,6\z)  are  normalized  such  that 


poo 

/  x¥m(z)x¥„(z)/p (z)dz  =  8mn 
Jo 


(13) 


and  assumed  to  vary  slowly  over  the  horizontal  plane  (r,  0);  but  the  complex  mode  amplitudes  Am{r,  Q) 
vary  much  more  rapidly  due  to  fluctuations  in  the  environment.  Following  Pierce  [32],  mode  coupling  is 
neglected. 

With  the  separation  of  variables  Eq.  (12),  the  Helmholtz  equation  Eq.  (11)  splits  into  a  Sturm-Liouville 
equation  Eq.  (14)  for  the  eigenfunctions  and  the  horizontal  wavenumbers  km(r ,  0);  and  a  two-dimensional 
Helmholtz  equation  Eq.  (15)  for  the  mode  coefficients  Am(r,6)  in  which  the  horizontal  wavenumbers  are 
effectively  the  index  of  refraction, 


P (z)  J;  ( d^m^dz  6 ' ^  ^  8 ; ' z)  ~  *"» ' (r, 6 ) )  (r,  d ; . z)  =  0,  and  (14) 


1  d  (  dAm(r,d)\  d2Am(r,6)  2  S(r) 


r  dr 


+  • 


,.2^02  +km(P°)Am(pO)  =  ~  Qs  m(0 , 0;  Zs)  ■ 


(15) 
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The  mode  coefficients  Am(r,d)  are  subject  to  stochastic  perturbations  due  to  sound-speed  inhomogeneities 
and  bottom  irregularities.  This  follows  from  k(r.  6;z)  =  G)/c(r,  9;z),  which  makes  km{r,  6 )  stochastic  in 
Eq.  (14),  and  Am(r ,  6)  stochastic  in  Eq.  (15).  Unlike  the  eigenfunctions,  these  environmental  perturbations 
are  assumed  to  vary  rapidly  over  the  horizontal  plane.  Because  normal-mode  amplitudes  propagate  in  the 
horizontal  plane  along  independent  paths  through  the  ocean,  they  are  affected  by  unrelated  sound-speed 
fluctuations.  This  process  causes  statistical  independence  among  the  mode  amplitudes. 

For  a  horizontal  line  array  oriented  broadside  toward  an  acoustic  source,  and  with  an  aperture  small 
compared  to  the  source-receiver  range,  curvature  of  the  acoustic  field  may  be  neglected,  and  the  receiver 
range  r  =  x„  and  depth  z0  are  constant.  These  assumptions  allow  a  simplification  of  notation:  y  =  xa9. 
P(y )  =  P{x0,y,Zo),  Am(y)  =  Am(x0,y),  and  'Fm(y)  =  lF„, (x0,y,z0)-  Furthermore  under  these  assumptions, 
eigenfunctions  that  vary  slowly  over  the  horizontal  plane  effectively  remain  constant  across  the  array  aper¬ 
ture,  vPm(y)  =  'Em,  My.  The  normal-mode  expansion,  Eq.  (12),  for  the  acoustic  pressure  field  evaluated  along 
a  broadside  horizontal  array  then  becomes 


M 


P(y)  =  Y  Am{yy¥m- 

m—  1 


(16) 


2.3  Modal  Expansion  of  Coherence  Function 

The  normal-mode  expansion  of  the  horizontal  coherence  function  follows  from  substituting  Eq.  (16)  into 
Eq.  (7), 


Rp(Ay) 


1  M  M 

-2  Y  I  <Am(y)A,Xv  +  Ay))t^n 

P  m—  1  n=  1 


(17) 


where  oj,  =  (|F(y)|2)fv  is  the  expected  value  of  the  acoustic  field  intensity.  The  mode  amplitudes  prop¬ 
agating  along  independent  paths  are  affected  by  unrelated  sound-speed  fluctuations,  hence  statistical  inde¬ 
pendence  of  the  fluctuations  make  the  cross-modal  amplitude  matrix  diagonal,  (Am(y)A*(y  +  Ay))t  ~  0, 
m  yk  n.  The  modal  expansion  of  the  cross-spectral  density  matrix  Eq.  (17)  then  becomes  a  single  sum  over 
the  modes, 


Rp(Ay) 


i  M 

-y  I  (A,n(y)K(y+ Ay))t,Xr 

°P  m=  1 


(18) 


Because  the  full  acoustic  field  is  homogeneous,  each  cross-modal  amplitude  matrix  is  also  homogeneous 
with  the  same  form  as  Eq.  (7), 


RAm(Ay)  = 


(Am(y)A*m(y  +  Ay))tiy 


(19) 


where  a\m  =  (Am  {y)2)t  is  the  expected  intensity  of  each  normal  mode,  and  R,\m  (Ay)  is  the  mutual-coherence 
function  of  each  mode  amplitude. 
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As  with  the  total  acoustic  field,  the  modal  coherence  functions  may  be  separated  into  an  amplitude 
correlation  and  a  modal  phase-structure  function  Dm  (Ay)  with  the  associated  normalized  form 

^Am(  Ay)  =  (Am(y)A*m(y  +  Ay))fy 

=  exp  [— ±Dm(Ay)]  ,  (20) 

where  Am(y)  =  Am(y)/|Am(y)|.  The  theory  developed  below  relates  Dm(Ay)  to  the  internal-wave  spectrum. 
The  normal-mode  expansion  of  the  horizontal  coherence  function  is  then 

M  sjl  m2 

RP( Ay)=t  ^V^A„,(Ay),  (21) 

m=  1  °P 

where  cry,  =  Y*m=\  shows  the  field  power  is  the  total  power  in  the  modes  evaluated  at  the  array  depth. 

Here  the  total  power  is  a  normalization  constant  that  enforces  Rp(0)  =  1. 

The  functional  form  of  Eq.  (21)  is  generally  true  for  any  statistical  process  that  is  a  sum  of  independent 
zero-mean  subprocesses  each  with  variance  a2.  A  similar  result  was  found  by  Shifrin  [35]  using  very 
general  statistical  arguments  for  radar-antenna  field  coherence.  Similarly,  for  an  acoustic  field  that  is  more 
appropriately  represented  by  a  supeiposition  of  rays,  the  horizontal  coherence  comprises  a  weighted  sum  of 
functions,  each  belonging  to  a  ray  that  connects  the  source  and  receiver. 

In  the  following  sections  a  path-integral  approach  is  used  in  conjunction  with  the  Markov  approximation 
and  a  spatial  spectrum  for  linear  internal  waves  to  obtain  an  expression  for  R,\m  (Ay). 


3,  MUTUAL  COHERENCE  OF  MODE  AMPLITUDES 

3.1  Path  Integral  Representation  of  Modal  Coherence 

To  obtain  a  path-integral  expression  for  the  horizontal  coherence,  the  wave  equation,  Eq.  (15),  for  the 
mode  amplitudes  is  converted  to  Cartesian  form.  Let  x  =  r,  the  cross-range  position  be  y  =  rd ,  and  drop  the 
explicit  dependence  on  horizontal  position,  (r,  6)  — >  ( x,y ),  from  the  notation  for  Am.  In  a  far-field  limit  the 
transverse  array  aperture  is  much  shorter  than  the  range,  i.e.,  y  <C  x0,  and  6  is  small  enough  that  wavefront 
curvature  may  be  neglected. 

Next  a  coordinate  transformation  is  introduced  that  renders  the  problem  Cartesian: 

Am  =  r-'/2Am,  (22) 

=  4 +  1/(2  rf.  (23) 

Henceforth  the  term  horizontal  wavenumber  refers  to  the  reduced  wavenumber  Km  defined  by  Eq.  (23). 

Clearly  this  approximation  breaks  down  where  r  <  1/2 km.  For  the  lowest  acoustic  modes  associated 
with  large  km,  breakdown  occurs  near  the  source.  The  path-integral  solution  presented  below  is  then  valid 
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from  a  few  wavelengths  beyond  the  source,  through  the  ocean  fluctuations  to  the  receiver.  For  the  highest 
modes  that  are  near  cutoff  where  km  —y  0,  the  path-integral  solution  breaks  down  over  the  entire  range. 

Thus  km  ~  Km  in  the  remainder  of  this  work.  With  this  change  of  spatial  coordinates  and  dependent 
variables,  it  can  be  shown  that  the  homogeneous  part  of  Eq.  (15)  reduces  to  the  Cartesian  form 


k)2A  ?2A 
~dxr  +  ~dy2~ 


+  ^i(x,y)Am 


=  0. 


(24) 


Following  Flatte  [10]  and  Dashen  et  al.  [27],  let  the  perturbations  of  the  horizontal  wavenumber  have  the 
form 

Kn(x,y)  =  km(l+um(x,y)  +  lJ,m(x,y)) ,  (25) 


where  km  =  ( Km(x,y))xy  is  the  mean  horizontal  wavenumber  over  the  entire  region  between  the  source  and 
receiver;  um(x,y)  accommodates  the  coordinate  transformation  and  describes  a  wavenumber  “profile”  that 
varies  slowly  over  the  horizontal  plane,  e.g.,  a  gradient  induced  by  propagation  on  a  slope;  and  fj,m(x,y) 
are  the  random  fluctuations  that  cause  statistical  variability  in  the  acoustic  field.  km  is  a  dimensional  scale 
factor  and  both  um  and  p,„  are  nondimensional.  An  important  assumption  here  is  that  both  the  deterministic 
perturbation  um  -C  1  and  stochastic  fluctuation  jlm  -C  1  are  small. 

For  the  propagation  of  mode  amplitudes  in  the  horizontal  plane,  the  dominant  ray  is  nearly  a  straight 
path  from  source  to  receiver.  Unlike  propagation  problems  in  the  vertical  plane,  there  are  no  multiple  paths 
here  except  near  caustics  that  form  at  mode  cutoff  for  upslope  configurations.  Thus,  it  is  reasonable  to  look 
for  solutions  of  the  form 

Am(x,y)  =  ym(x,y) exp  ( ikmx ) .  (26) 


In  the  small-angle  parabolic  approximation  of  Eq.  (24),  the  d2\ \rm/dx2  ~  0  term  is  ignored,  as  are  the 
quadratic  wavenumber  perturbations  u2n,  /A2,  and  thus  producing 


2 1  Kjfi 


d  yu 

dx 


d2Wm 
dy 2 


2tcm  (um (x,y)  -F l~lin (x.y) )  \(/m . 


(27) 


(Also,  a  pair  of  low-order  k2  \ jfm  terms  of  opposite  sign  cancel.)  As  Tatarskii  [36,  p.  375]  notes,  the  parabolic 
approximation  is  valid  for  small  wavelengths  where  both  A,„  -C  le  and  X2nx0  <C  /f4.  Here  X,„  =  2n/km  is  the 
horizontal  wavelength  of  each  mode;  le  is  the  scale  length  of  environmental  fluctuations;  and  xQ  is  the 
propagation  range.  These  criteria  clearly  break  down  as  km  — >  0,  which  occurs  for  high-order  modes  near 
cutoff. 


The  path-integral  solution  to  Eq.  (27)  is  given  by  Flatte  [10],  Dashen  [11],  Esswein  and  Flatte  [12],  Flatte 
et  al.  [25],  and  Dashen  et  al.  [27]: 


1/2 


Dp  exp  <  iK, 


[Um  (x,  yp  (x) )  +  ixm  (x,  yp  (x) )]  dr 


(Xq :  Y )  —  (l/4Km) 


(28) 
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Um(x,yp(x)) 


1  /  dyp(x) 

2  y  dx 


2 

+  um{x,yp{x)), 


(29) 


where  the  subscript  p  denotes  the  path  and  fDp  indicates  integration  over  all  possible  paths.  In  the  above 
integral,  all  paths  begin  at  the  origin  and  terminate  at  the  same  point,  Y  =  yp{x0),  V/;.  Here  Um(x.yp(x))  pro¬ 
vides  the  solution  in  the  absence  of  environmental  fluctuations,  while  fxm(x,yp(x))  introduces  contributions 
from  the  fluctuations. 


Dashen  [11]  and  Dashen  et  al.  [27]  have  demonstrated  that  the  second  moment  of  Eq.  (28)  may  be  written 
as  an  exponential  of  the  phase- structure  function, 


(Vm(.Xo,YaWm(Xo,Yp))e  {  ,  , 

^ - =  exp  [-%Dm(Ya,Yp)\ 


{Xoi  0)  | 

where  the  path-integral  solution  for  the  structure  function  is 


(30) 


Dm(Ya,Yp)  =  k-{  (  j  Hm(xi,ya{xi))dxi- J  ) ■im(x2,yp(x2))dx2 


(31) 


The  points  Ya  =  ya(x0)  and  Yp  =  yp(xa)  are  the  terminal  transverse  coordinates  of  two  nearby  principal 
rays  a  and  /J  propagating  through  the  unperturbed  environment.  The  geometry  of  the  principal  rays  and 
the  coordinate  system  describing  them  is  shown  in  Fig.  1.  The  expectation  is  taken  over  an  ensemble  of 
homogeneous  stochastic  environments  e,  thus  (J  \jfm (x,y)|2^)  =  ^ |  !//,„ | 2 ^  for  all  x.y.  Also,  use  was  made  of 

the  identity  (exp(ia))  =  exp(— ^  (a2)),  which  is  valid  for  Gaussian  random  variables  that  result  when  the 
central-limit  theorem  is  applied  to  a  =  f  pm(x)dx,  even  for  non-Gaussian  pm(x). 


It  must  be  emphasized  here  that  Eq.  (31)  is  derived  from  a  double  path  integral  obtained  by  integrat¬ 
ing  solutions  to  the  parabolic  equation  over  all  possible  paths  between  a  source  and  two  receivers  ( a  and 
/3).  These  paths  through  a  stochastic  environment  fluctuate  about  two  deterministic  principal  rays  from  the 
source  to  the  pair  of  receivers.  The  principal  rays  are  defined  within  the  unperturbed  environment,  and 
diverge  approximately  linearly  with  range  as  discussed  in  Sec.  3.3.  From  the  integrals  of  the  index  of  re¬ 
fraction  over  range  in  Eq.  (31),  the  phase  structure  is  seen  to  arise  from  the  expected  value  of  stochastic 
differences  in  path  length,  and  has  the  same  form  as  Eq.  (3)  with  a  constant  wavenumber. 


These  results  for  mutual  coherence  are  valid  regardless  of  scattering  strength  up  until  the  parabolic  ap¬ 
proximation  breaks  down,  and  are  not  based  on  a  stationary  phase  approximation  through  the  fluctuations 
[27].  In  the  path-integral  formulation  of  mutual  coherence,  global  statistical  properties  reduce  the  integration 
over  all  possible  paths  to  an  integration  over  a  single  mean  path  [11].  The  conditions  of  validity  for  Eqs.  (30) 
and  (31)  are  rigorously  demonstrated  in  Appendix  A  of  Ref.  11.  The  main  result  is  that  for  an  isotropic  ho¬ 
mogeneous  medium,  Eq.  (30)  is  valid  when  the  parabolic  approximation  is  valid.  The  path-integral  solution 
is  very  general  [11,  27],  The  assumptions  are  the  parabolic  approximation,  homogeneity  of  environmental 
fluctuations  in  the  transverse  direction,  and  the  existence  of  a  single  principal  ray  arriving  at  each  receiver. 
Thus,  these  results  are  not  valid  in  a  multipath  environment,  such  as  shallow-water  acoustic  propagation  in 
a  vertical  plane. 
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Fig.  1  —  Geometry  of  two  principal  rays  a  and  j8  propagating  in  the  horizontal  plane  (x,y).  The 
acoustic  source  is  located  at  the  origin,  and  the  transverse  receiver  array  represented  by  the  vertical 
line  of  dots  is  located  at  range  x0 .  The  nested  ellipses  represent  the  two-point  correlation  function, 
pm(Ax,Aap)  =  ( jlm (x\  ,y«) Pm {x2 .  >yi ) ) ,  °f  fluctuations  in  the  horizontal  acoustic  wavenumber 
A lm(x,y). 


For  our  problem  of  mode  amplitudes  propagating  in  a  horizontal  plane  without  lateral  boundaries,  the 
assumptions  required  by  path-integral  theory  arc  reasonable.  Here  there  are  no  lateral  boundaries  to  cause 
multiple  paths.  Mode  amplitudes  propagating  in  the  horizontal  plane  governed  by  Eq.  (24)  follow  benign 
paths  that  diverge  slowly,  except  up  a  slope  near  mode  cutoff  where  caustics  form  [33,  37,  38].  It  remains 
now  to  determine  precisely  the  decay  rate  for  homogeneous  isotropic  shallow-water  environments. 

3.2  Markov  Approximation  of  Phase-Structure  Function 

Following  Esswein  and  Flatte  [12],  the  phase-structure  function  Eq.  (31)  is  expressed  using  a  weak  form 
of  the  Markov  approximation.  Unlike  stratified  propagation  problems  formulated  in  the  vertical  plane,  here 
there  is  no  anisotropy  to  contend  with,  which  simplifies  the  final  results. 

Equation  (31)  is  insensitive  to  interchanging  the  two  paths,  and  they  may  be  combined  when  the  square 
is  expanded, 

rxo  rxo 

Dm(Ya,Yp)  =  2K%lJ  Pm{xi,x2;ya(xi),ya(x2))-Pm{xi,x2;ya(xi),yp{x2))dxidx2,  (32) 
where  the  two-point  correlation  function  of  the  environment  between  principal  rays  a  and  /J  is 


Pm  {x\ ,  x2;ya  (xi ) ,  yp  {x2))  =  (pm  (u ,  ya  (*i ) )  Pm  (x2,yp(x2)))e- 


(33) 
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Transverse  homogeneity  was  assumed  in  the  derivation  of  Eqs.  (30)  and  (31),  and  here  homogeneity  in  range 
is  assumed  as  well: 

Pm(xi,X2;ya(xi),yp(x2))  =  Pm(x2-xi,yp(x2)  -ya(xi))-  (34) 

Additionally,  the  range  variables  are  changed  to  a  range  difference  and  a  mean  range:  Ax  =  x2  —  x\,  x  = 
{x\  +V2)/2;  with  associated  changes  in  the  cross-range  direction:  Av  =  Aap(x0),  Aap (x)  =  yp (x)  ~~ya[x). 
Again,  see  Fig.  1  for  a  sketch  of  the  geometry. 

Now  the  Markov  approximation  is  made.  Because  pm(Ax,Aap(x))  is  compact  around  Ax  =  0,  the  corre¬ 
lation  function  has  support  in  range  only  where  x\  ~  X2  ~  x,  and  is  negligible  elsewhere.  Hence,  Esswein 
and  Flatte  [12]  justify  the  approximation 

Pm(x2- xi, yp(x2)  -ya(xi))  & pm( Ax,Aap(x)).  (35) 

Compactness  also  justifies  extension  of  the  integration  limits, 

r*o  r°° 

j  Pm(Ax,Aap  (x))dAx  J  pm (Ax.:  Aap  (x j  jdA, .  (36) 

In  the  strong  Markov  approximation  the  correlation  in  range  is  infinitesimal, 

Pm{Ax,Aap(x))  ~  8 {Axs)pm(^Aap (v)).  (37) 

With  the  weak  Markov  approximation,  the  correlation  in  range  remains  finite.  With  these  homogeneity  and 
compactness  assumptions,  the  phase-structure  function,  Eq.  (32),  becomes 

An  (Ay)  =2Kmfo  J  Pm(Ax,0)-pm(Ax,Aap(x))dAxdx.  (38) 

The  Markov  approximation  breaks  down  near  turning  points  where  the  ray  curvature  is  large  in  com¬ 
parison  with  the  correlation  length  of  environmental  fluctuations.  This  difficulty  does  not  occur  for  the 
shallow-water  propagation  of  normal-mode  amplitudes  in  the  horizontal  plane.  Here,  the  problem  of  large 
ray  curvature  occurs  only  on  slopes  in  the  vicinity  of  caustics  near  mode  cutoff. 

3.3  Spectral  Representation  of  Phase-Structure  Function 

The  stochastic  acoustic  horizontal-wavenumber  fluctuations  pm(x,y)  are  zero  mean  by  definition,  and 
in  the  development  of  the  phase-structure  function  the  process  was  also  assumed  homogeneous.  It  is  thus 
reasonable  to  assume  that  the  fluctuations  belong  to  a  stationary  stochastic  process.  Under  these  conditions 
the  correlation  function  pm(Ax.Aap(x))  may  be  represented  in  spectral  form  via  the  Wiener-Khintchine 
theorem  [39], 

pm(Ax,Aap(x))=  [  [  Sm(yx,Yy)eil(7A+^AaP{-x))Uyxdyy,  (39) 

where  yx  and  yy  are  the  horizontal  wavenumbers  in  the  range  and  transverse  directions  for  the  power  spectrum 
Sm(Yx,Yy )  °f  the  acoustic  horizontal-wavenumber  fluctuations.  Note  there  are  two  very  different  types  of 
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horizontal  wavenumber  in  this  problem.  The  horizontal  wavenumbers  of  the  acoustic  field,  Km,  act  as  indices 
of  refraction  governing  the  horizontal  propagation  of  the  normal-mode  amplitudes,  whether  the  propagation 
be  deterministic  or  stochastic.  When  the  propagation  is  stochastic,  fluctuations  of  this  index,  pm,  also  have 
a  power  spectrum  S  with  horizontal  wavenumbers  denoted  by  yx  and  yv. 

The  integrals  that  compose  the  phase-structure  function  Eq.  (38)  are  represented  in  spectral  form  as 

/oo  poo  poo  .  . 

pm(Ax,Aap(x))dAx  =  2n  /  Sm(yx,yy)  8(yx)eV7yAa^x))dyxdyy,  (40) 

-oo  J — oo  J — oo 


after  substituting  dAx  =  2n8(yx).  Now  assume  horizontal  isotropy  by  letting  yx  =  y,  cos(0)  and 

Yy  =  Yu sin(0),  and  note  that  with  this  change  of  coordinates  5(yv)  — >  8(y/,cos(0j)  and  dyvdyv  H>  y,dy,dd 
where  the  delta  function  is  nonzero  for  y,  >  0  when  6  =  2 mi  ±n/2,  then  the  correlation  function  becomes 


Pm  (AXlAap  (x)  )dAv  —  In 


e{+iyiAap{xj)  _|_g(-'r;,Aa/5W) 


Yhdyh, 


(41) 


and  likewise  when  =  0, 


pm(Ax,0)dAx  =  An 


Sm(y,)yAYh- 


(42) 


Because  of  the  properties  of  8(yl  cos(0)),  yx  =  0  and  yy  =  ±y„  thus  showing  that  Eqs.  (41)  and  (42)  depend 
solely  on  the  transverse  part  of  the  spectrum.  From  this  it  is  apparent  that  only  environmental  fluctuations 
with  wavevectors  normal  to  the  main  acoustic  path  contribute  to  the  transverse  coherence.  Contributions 
from  fluctuations  with  wavevectors  parallel  to  the  acoustic  path  average  out  to  zero.  The  spectral  represen¬ 
tation  of  the  phase-structure  function,  Eq.  (38),  then  becomes 


DjAy)  =4KKfn  /  /  sm(yh) 

Jo  Jo 


2- 


?(+iVlAall(x))  _|_  e(-iYhAap(x)) 


y,dy,dx. 


(43) 


To  evaluate  the  integral  over  range  x,  make  the  approximation  that  two  nearby  principal  rays  a  and  /3  di¬ 
verge  approximately  linearly,  Aap  (x)  ~  Ayx/x0,  as  they  propagate  in  the  horizontal  plane  of  the  unperturbed 
environment.  This  approximation  becomes  exact  for  a  constant-depth  environment  with  uniform  sound 
speed.  Even  when  the  rays  curve  weakly  on  a  slope  far  from  caustics,  their  relative  divergence  remains 
approximately  linear  in  range.  With  this  approximation  the  integral  may  be  evaluated, 


eO'vA  a/>W)dx=(x0/iy,Ay) 


e(‘My)  _  l 


(44) 


and  the  phase-structure  function  becomes 


JJm  (Ay) 


poo 

4kx0kI  /  Sm(yh) 
Jo 

e(ij,Ay)  _  e(-iTlAy) 

YhdYh 

(45) 

iYh^y 

poo 

8ux0kI  /  Sm(yh)  [1  -  sine  (%Ay)]  yhdyh. 
Jo 

(46) 
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It  must  be  emphasized  here  that  the  sine  function  arises  from  an  integral  over  range  of  phase-difference 
fluctuations  between  two  principal  rays,  which  happen  to  diverge  linearly.  It  has  nothing  to  do  with  restrict¬ 
ing  integration  over  all  possible  paths  to  a  finite  transverse  region.  Integration  over  all  possible  paths  was 
done  properly  in  Sec.  3.1. 

The  unknown  spectral  weighting  Sm{%)  makes  it  impossible  to  say  much  about  the  integral  except  for 
limiting  cases.  A  power-series  expansion  of  the  sine  function  for  small  Av  has  quadratic  behavior, 

1  -  sine (y,, A y)  «  (y/,Av)2/6  -  <^(y/,Ay)4.  (47) 

Conversely  in  the  opposite  limit  Ay  — >•  °o,  the  sine  function  decays  as  the  phase- structure  function  approaches 
the  spectral  form  of  the  modal-scattering-strength  parameter 

poo 

Dm{oo)  =  %TzXokl  Sm(jh)  YiAYh,  (48) 

7o 

where  the  integral  is  the  total  power  in  the  horizontal-wavenumber  spectrum. 

The  spectrum  Sm(y/,)  required  by  the  phase- structure  function  can  be  obtained  from  a  known  correlation 
function  pm{Ax.Ay).  Because  the  internal-wave  displacement  fluctuations  are  assumed  homogeneous  and 
isotropic,  the  correlation  function  may  be  evaluated  at  any  point  in  the  ocean  along  any  orientation.  Thus  it 
is  sufficient  to  evaluate  Eq.  (41)  at  x  =  xa,  where  Aap  (x)  =  Ay,  and  solve  for  STO(%).  Equation  (41)  may  then 
be  rewritten  as 

/oo  poo 

pm(Ax,Ay)dAx  =  2n  Sm{yh)eiyh^\'Yh\d'Yh,  (49) 

-oo  J  — oo 


where  symmetry  of  the  power  spectrum  Sm(jh)  has  been  exploited  to  combine  the  two  halves  of  the  integral. 
Equation  (49)  is  then  Fourier  transformed  to  obtain 

1  poo  poo 

Sjjh)  =  2n\n\  J  J  P,n(Ax.Ay)e-iV‘A>dAxdAy.  (50) 

To  make  progress,  either  the  spectrum  or  the  correlation  function  must  be  connected  to  the  underlying 
environmental  fluctuations  by  direct  experimental  measurement,  or  by  theoretical  analysis  of  an  internal- 
wave  model.  The  latter  approach  is  pursued  in  the  following  sections. 

The  appropriate  fluctuation  spectrum  for  shallow  oceans  is  neither  well  established  nor  universal.  Con¬ 
tributory  processes  that  induce  randomness  include  the  thermohaline  microstructure,  linear  and  nonlinear 
internal  waves,  bottom  impedance,  bottom  roughness,  and  surface  waves  [40].  These  processes  are  site  de¬ 
pendent  and,  excepting  bottom  properties,  influenced  by  the  seasonal  stratification,  tidal  cycle,  and  weather. 
To  proceed,  the  theoretical  analysis  of  horizontal  coherence  must  be  restricted  to  a  specific  environment. 
Among  the  above  processes,  shallow-water  environments  often  show  evidence  for  the  impact  of  internal 
waves,  the  process  that  will  be  pursued  here. 
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4.  ACOUSTIC  AND  INTERNAL- WAVE  SPECTRA 

In  this  section,  a  relationship  between  fluctuations  in  the  acoustic  horizontal  wavenumbers  and  those 
of  the  internal-wave  spectrum  is  obtained  based  on  a  perturbation  expansion  of  the  Helmholtz  equation,  in 
which  the  wavenumbers  are  linearized  for  small  perturbations  of  the  sound  speed  field.  The  sound-speed 
fluctuations  themselves  are  related  to  internal-wave  displacements  using  the  potential  sound-speed  gradient 
and  linear  internal-wave  modes  as  basis  functions.  Correlation  functions  are  formed  by  taking  the  expected 
value  of  fluctuation  products,  and  Eq.  (50)  is  used  to  obtain  the  spectral  relationship.  Finally,  a  Garrett-Munk 
spectrum  modified  for  shallow  water  is  used  for  the  transverse  spatial  spectrum  of  the  internal  waves. 

4.1  Sound-Speed  Perturbations 

In  App.  A,  a  general  solution  is  obtained  for  the  relationship  between  stochastic  fluctuations  in  the  eigen¬ 
values  of  the  Helmholtz  equation  and  those  of  the  sound-speed  field.  For  the  acoustics  problem  under 
consideration,  the  solution  for  the  eigenvalue  fluctuations  given  by  Eq.  (A  15)  has  the  form 

—  9 

Hm(x,y)  =  1_ J  v(x,y,z)^,n(z)/p{z)dz.  (51) 

Ki  Jo 


Equation  (51)  asserts  that  fluctuations  in  the  horizontal  wavenumber  arc  approximately  proportional  in  an 
integral  sense  to  fluctuations  in  the  total  acoustic  wavenumber.  The  stochastic  sound-speed  perturbations 
v(x,y;z)  required  by  Eq.  (51)  are  obtained  from  internal-wave  displacements.  In  Sec.  4.3  below,  Eq.  (51)  is 
squared  and  an  expected  value  taken  to  obtain  pm,  which  is  required  by  Eq.  (33)  for  use  in  Eq.  (50). 


4.2  Internal- Wave  Displacements 


In  a  stratified  ocean,  the  sound-speed  fluctuations  within  the  water  column  are  dominated  by  linear 
internal-wave  displacements  [9,  41].  Under  these  conditions  the  sound-speed  gradient  dc/dz  =  dpc/dz  + 
dac/dz  is  the  sum  of  an  adiabatic  gradient  da  and  a  potential  gradient  dp,  where  only  vertical  fluctuations 
within  the  potential  gradient  contribute  to  the  sound-speed  fluctuations  [8,  42,  43].  A  relationship  between 
vertical  displacements  and  sound-speed  fluctuations  is  then  obtained  by  equating  first-order  terms  from  the 
Taylor-series  expansion  of  c{z)  with  similar  terms  from  the  sound-speed  perturbation  expansion  c(z)  = 
c(l  —  v(z)  —  v(z))  to  obtain 


v(x,y;z) 


1  dpc 
c  dz 


(x,y,z)^(x,y\z). 


(52) 


Following  Flatte  and  Tappert  [8],  Gill  [44],  and  Tielbiirger  et  al.  [45],  the  vertical  displacements  of  the 
water  column,  £,  are  expanded  as  a  sum  of  J  linear  internal-wave  modes  Wj,  with  complex  mode  amplitudes 

8j> 


N  pco  p27 1 

CO x,yz,t)  =  J^  / 

Pi  Jo  Jo 


g jilt'll,  d)Wj(x,y;z,  7/i)e!^K!COs(0)+:yl,'sin(0)_COjf]d0  7/jdy/,. 


(53) 


where  (x,y)  are  horizontal  coordinates,  t  is  time,  and  z  is  depth  with  the  free  surface  at  z  =  0  and  ocean 
bottom  at  z  =  z,„-  The  associated  Cartesian  wavenumbers  arc  yx  =  y/,cos(0)  and  yv  =  y, sin(0),  which  in 
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cylindrical  coordinates  in  wavenumber  space  are  written  in  terms  of  the  total  horizontal  wavenumber  //,  and 
an  angle  in  the  horizontal  plane  6. 

The  vertical  oscillation  frequencies  (Oj(Yh)  and  internal-wave  mode  functions  Wj{z,%)  are  the  param¬ 
eterized  eigenvalues  (dispersion  relation)  and  eigenfunctions  of  a  well-known  Sturm-Liouville  problem 
[25,  44,  46-49]  for  the  unperturbed,  stratified,  spatially  homogeneous,  mean  environment.  The  Wj(z,  Yh)  are 
deterministic  mode  functions,  with  associated  buoyancy  frequency  N(z),  of  a  stratified  spatially  homoge¬ 
neous  mean  environment  that  satisfies 


d2Wj(z,Yh )  ,  ,2 

dz2  +Yfl 


N2(z )  -  co2 


Wj(z,Yk)=  0, 


(54) 


such  that  Wj(0)  =  0,  Wj(za)  =  0,  and  with  a  normalization  that  ignores  the  Earth’s  rotation  [25], 

N2(z)Wi(z,  Yh)Wj(z,  Yk)dz  =  8U,  (55) 

Jo 

where  is  the  Kronecker  delta.  Here  f,  =  2£2sin(0)  is  the  local  Coriolis  frequency,  Q.  =  7.292  10_5rad/s 
is  the  Earth’s  rotation  rate,  and  0  is  the  latitude. 

The  mode  coefficients  gj(Yh-O)  are  samples  of  stochastic,  zero-mean,  complex  amplitudes  [17,  25] 
drawn  from  a  distribution  having  a  second  moment  <  gjg*j  >e  that  is  isotropic  and  homogeneous.  It  is 
important  to  include  the  angular  dependence  in  gj(%,6)  because  internal  waves  arriving  from  different 
directions  are  generated  by  unrelated  physical  processes. 

Unless  the  bottom  slope  or  horizontal  density  gradient  is  large,  e.g.,  in  the  vicinity  of  a  shelf-break 
front,  the  unperturbed  seasonal  stratification  is  nearly  uniform  over  a  horizontal  plane  that  spans  just  a  few 
tens  of  kilometers.  Under  these  conditions  there  is  no  coupling  between  internal-wave  modes,  and  the 
associated  mode  functions  are  approximately  independent  of  horizontal  position,  Wj(x,y;z ,  Yh)  ~  7/,), 

N(x,y,z)  ~N(z),  and  dpc(x,y,z)/dz  ~  dpc(z)/dz. 

With  these  approximations,  Eqs.  (51) -(53)  become 


-k2  N 

i Um{x,y,t)  =  — z^-  52 
j= f 


r2n 


1C^ 


^mjijYh) 


gy(7/„  0)c'[xr',cos(0)+v7'lSin(e)”H de  YhdYh, 


and 


(56) 


\^{z)^m{z)wj{z,Yh)dz. 


(57) 


The  importance  of  Eq.  (57)  is  that  it  weights  most  heavily  those  internal-wave  displacements  at  depths 
where  the  potential  sound-speed  gradient  is  large,  i.e.,  within  the  thermocline,  and  deemphasizes  high- 
order  internal  modes  that  oscillate  rapidly  in  that  part  of  the  integral.  The  impact  of  this  process  on  any 
particular  acoustic  mode  depends  on  how  the  shape  of  the  acoustic  mode  function  aligns  with  the  internal- 
wave  excitation.  Acoustic  modes  with  maxima  in  the  vicinity  of  large  potential  sound-speed  gradients  suffer 
more  severe  transverse  coherence  loss  than  modes  with  a  null  near  the  gradient. 


16 


Peter  C.  Mignerey 


4.3  Internal- Wave  Correlation  Function 

A  relationship  between  the  acoustic  correlation  function  and  the  internal-wave  spectrum  is  obtained  by 
substituting  Eq.  (56)  into  Eq.  (33): 


T.2  \  2  N  N 
Pm(  Ay,Av)  =  (  —  )  £  £ 


poo  poo 


j  1=1  72=1 


1 0  JO 


-  Kri.j]  (Y/i i  )-Kn.ji  (.Yli2 


p2k  r 2k 


10  Jo 


(gjl(lkl,0i)g*j2(Yk1,02)) 


A*1* i  c°s(0i)-t-yi^i  sin cosCftJ+wftj  sin(e2)-0L.72r] d0i d02  ^  y^dy^dy,,, .  (58) 


The  statistical  properties  of  internal-wave  modes  are  now  exploited  with  respect  to  mode  number,  wavenum¬ 
ber,  and  propagation  direction.  First,  internal  waves  arriving  from  different  directions  are  generated  by 
unrelated  physical  processes  that  are  widely  separated  spatially.  Therefore  such  arrivals  are  assumed  to 
be  statistically  uncorrelated,  (g-n  (//,, ,  B\  )g*j2{Yh2 ■  B^'ij  =0,  B\  y  Bi.  Second,  second-order  statistics  of 

the  generating  mechanisms  are  assumed  to  be  uniform  in  azimuth,  ^  \gj(Yh-  0 )  |2^)  =  Gj{%),  V0.  Third,  be¬ 
cause  the  internal-wave  modes  are  uncoupled,  the  mode-amplitude  fluctuations  are  statistically  independent, 
(^gji  (yhl ,  0i  )g*j2  (jh2  ,  0t)^  =0,  yi  /  yV  Fourth,  linear  internal-wave  amplitude  fluctuations  are  assumed  to 

be  governed  by  an  orthogonal  stationary  stochastic  process  [39],  (gjiYh, , B\  )g*j{Yh2,  B^'i'j  =  0,  Yh  y  y/,2. 
With  these  statistical  assumptions,  the  internal-wave  power  spectrum  simplifies  to 

(8jSYhl,0i)g*2(Yh2,e2))e  =  Gj(Yh)8jlj28{Yhl-Yh2)8^eiv  &2\  (59) 

Yh 


and  the  correlation  function  of  the  horizontal  wavenumbers  becomes 


Pm  (Ax  1  Ay  ) 


eiyh  A-'  «*(»)  giftAr  si"( e)  d  0  y/(dy/( . 


(60) 


4.4  Internal- Wave  Spectrum 

The  relationship  between  acoustic  and  internal-wave  spectra  is  finally  obtained  by  substituting  Eq.  (60) 
into  Eq.  (50), 


Sm  ( Yh  ) 


piYh'^c  osWdA 


errh'Ay  M8)e-in A dAydB  y/dy,/ . 


(61) 


Note  here  that  Ju  >  0  because  of  the  limits  of  integration.  Thus  Yh'  may  be  replaced  by  Yh'  within  the 
integral.  Then  recognizing  that 


/oo 

govA,cos(e)dAc  =  27r5(yi,Cos(0)) 

-00 


(62) 
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and  that  5(y//cos(0))  is  nonzero  for  Yh'  >  0  only  when  6  =  2nn  +  n/2,  transverse  wavenumbers  are  again 
selected;  likewise: 


/oo 

e±,ri''Aye-,Yl'A>'dAy  =  2nS(±yy  —  //,). 

-oo 


(63) 


With  these  observations,  the  spectral  relationship  Eq.  (61)  becomes 


2n  (  k2  \  2  N  f°° 

Sm(Yh)  =  7771(^2)  L  /  A nj(Yh')Gj{Yh >)  [8{Yh’  -  Yh)  +  8{yw  +  Yi ,)]  IM&Yh'-  (64) 

\Yh\  \  /  j=\  70 


The  properties  of  the  delta  functions  imply  //,  =  ±y/,/.  But  because  //,,  >  0,  one  delta  function  defines 
the  positive  half  spectrum  Sm(Yh  >  0),  and  the  other  the  negative  half  spectrum  Sm(Yi,  <  0).  Although  the 
Sm(Yh  =  0)  component  is  counted  twice,  it  is  zero.  Hence  the  spectrum  simplifies  to 

Sm{Yh)  =  ^n(y\  Z<j(n)Gj{n).  (65) 

V  Km  J  j—  1 

Then  substituting  Eq.  (65)  into  Eq.  (46),  an  expression  for  the  modal  phase- structure  function  is  obtained  in 
terms  of  the  internal-wave  spectra, 


^m,j{Yh)Gj{Yh)  [1 


sine  (Yh  Ay)]  y/,dy,,, 


(66) 


where  Hm  =  1 6n2x(lk4  /  Kyn. 

An  observation  of  the  coherence  dependence  on  acoustic  mode  number  may  now  be  made.  Because 
the  horizontal  acoustic  wavenumbers  km  decrease  with  increasing  acoustic  mode  number  m,  the  phase- 
structure  function,  Eq.  (66),  increases  with  mode  number  causing  the  exponential  in  the  modal  coherence 
function,  Eq.  (20),  to  fall  off  faster  with  transverse  separation  for  the  higher  modes  than  for  the  lower  modes. 
Likewise,  Hm  increases  with  total  wavenumber  k  as  the  frequency  rises,  again  causing  a  loss  of  coherence. 
It  remains  now  to  introduce  an  appropriate  shallow-water  model  for  the  internal-wave  spectrum  G;(y/,). 

4.5  Internal- Wave  Spectral  Model 

A  definitive  description  of  the  internal-wave  spectrum  in  shallow  water  remains  an  open  question.  Thus 
any  choice  is  somewhat  arbitrary.  One  candidate  is  a  modification  of  the  Garrett-Munk  spectrum  [50]  for 
shallow  water.  While  the  mechanism  for  the  generation  of  nonlinear  internal  waves  by  tidal  flow  across 
the  continental  shelf  break  is  well  established  [22,  51,  52],  the  generation  of  diffuse  shallow- water  internal 
waves  is  less  well  understood.  One  source  of  such  waves  is  the  arrival  from  afar  of  diffuse  deep-water 
linear  internal  waves  that  were  generated  in  the  open  ocean.  This  latter  process  provides  some  justification 
for  adapting  the  Garrett-Munk  spectrum  to  shallow  water  as  was  done  by  Yang  and  Yoo  [53].  However, 
it  must  be  kept  in  mind  that  both  the  deep-water  Garrett-Munk  spectrum  and  especially  any  shallow-water 
modifications  are  based  on  empirical  analyses. 
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Following  the  standard  definition  of  the  Garrett-Munk  spectrum  [41,  50],  and  using  Eq.  (53)  with 
Eqs.  (59)  and  (55)  to  obtain  the  variance  of  internal-wave  displacements  (J £|2^,  the  energy  per  unit  area 
of  ocean  floor  $  (J/ m2 )  is  given  by 

&  =  P  Jq  N2(z)(\£(x,y;z,t)\2^dz  (67) 

N  poo 

=  2 7Tp£  /  Gj(yh)  Yhdyh,  (68) 

H i7° 

where  p  is  a  global  mean  density  scale.  In  the  same  references,  the  Garrett-Munk  spectrum  £)(//,)  is  defined 
to  satisfy 

N  poo 

S=YJ  Ej(Yh)dYh.  (69) 

(Note:  The  Garrett  and  Munk  spectrum  in  Flatte  et  al.  [25,  p.  56]  includes  the  density  in  their  Eq.  3.2.22  for 
the  internal-wave  eigenfunction  normalization,  whereas  here  the  standard  normalization  given  by  Eq.  (55) 
does  not.)  Thus  the  relationship  between  the  Garrett-Munk  spectrum  and  the  internal-wave  spectrum  as 
defined  in  this  paper  is 

Ej(Yk)=2npYhGj(Yh)-  (70) 

The  Garrett-Munk  spectrum  of  each  internal-wave  mode  has  the  form 

Ej(Yh)=E0MjEj{Yh),  (71) 

where  Eq  is  the  total  internal-wave  energy  density  in  shallow  water,  EqMj  is  the  energy  of  each  mode,  and 

Mj  =  M(j2  +  jl)~p/2  (72) 

is  the  normalized  partition  of  energy  among  modes  with  constant  M  defined  such  that  Ylj=\Mj  =  1.  Rea¬ 
sonable  estimates  of  the  constants  are  j *  =  1  and  p  =  4  for  shallow-water  internal  waves  [53].  For  the 
deep-water  Garrett-Munk  spectrum  [9],  the  constants  are  j*  =  3  and  p  =  2. 

EjiYh)  is  either  the  normalized  deep-water  Garrett-Munk  spectrum 


EjiYh)  — 


or  the  Yang-Yoo  version  modified  for  shallow  water  [53], 


EjiYh)  = 


QjYjYl 


(jl  +  iYjN/fc)2)  V2 


(73) 


(74) 
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The  normalization  constant  Qj  is  chosen  such  that  J(“  Ej(jh)^7h  =  1-  For  the  shallow-water  form,  fc  is 
the  Coriolis  frequency,  the  wavenumber  spectrum  peaks  near  Jj  =  njfc/Nzo,  and  the  average  buoyancy 
frequency  is  defined  as  N  =  zG  1  f0"  N (z)d~. 

Substitution  of  the  Garrett-Munk  internal-wave  spectrum,  Eqs.  (70)  and  (71),  into  Eq.  (65)  produces  a 
Garrett-Munk  acoustic-mode  fluctuation  spectrum, 


5/;/  (  Yli ) 


Eo_ 

PYh 


Mj^nj(rh)Ej{yh). 


(75) 


Likewise,  the  phase-structure  function,  Eq.  (66),  in  either  deep  or  shallow  water  has  the  form 

Dm(Ay)  =  £  MjSmJ(Ay),  with  (76) 

^ P  j=  1 

ooo 

jYmJ(Ay)  =  jV^j{Yh)E ’j(jh)  [1  —  sinc(y/,Av)]  d//,.  (77) 

Immediately  apparent  is  the  direct  dependence  of  the  phase- structure  function  on  Eq,  which  causes  an  as¬ 
sociated  loss  of  transverse  coherence  in  the  presence  of  energetic  internal  waves.  The  dependencies  on  Mj 
and  -yVm.j  imply  the  impact  on  acoustic  coherence  is  most  severe  for  the  lowest  internal-wave  mode.  It  is  the 
integral  (Av)  that  determines  the  shape  of  the  phase  structure  as  a  function  of  Av  through  the  behavior 
of  [1  —  sinc(y/,Av)]  as  weighted  by  the  spectral  properties  of  the  internal  waves  and  buoyancy. 


5,  SIMULATED  COHERENCE 

To  illustrate  the  behavior  of  the  path-integral  solution  for  the  phase-structure  function  and  the  associated 
acoustic-field  coherence  that  was  developed  in  Ch.  3,  it  is  necessary  to  obtain  a  somewhat  realistic  power 
spectrum  for  the  eigenvalue  fluctuations.  The  required  spectrum  is  obtained  by  finding  parameters  for  the 
Garrett-Munk  spectrum  of  Sec.  4.5  that  produce  a  reasonable  fit  to  measured  shallow-water  internal-wave 
spectra. 

5.1  Experimental  Spectra 

For  the  above  purposes,  three  spectra  were  used  that  were  obtained  by  experiments  in  the  Barents  Sea 
(Barents)  [54,  55],  the  Mid-Atlantic  Bight  (SWARM)  [23],  and  the  East  China  Sea  (TAVEX)  [56,  57].  These 
environments  provided  a  diverse  set  of  spectra,  with  Barents  emphasizing  low  wavenumbers,  SWARM  mid 
wavenumbers,  and  TAVEX  high  wavenumbers. 

5. 1. 1  SWARM  and  Barents  Spectra 

Some  preliminary  comments  concerning  the  SWARM  spectra  are  necessary.  The  SWARM  spectra  were 
estimated  from  vertical- velocity  fluctuations  obtained  from  a  bottom-mounted  acoustic  Doppler  current  pro¬ 
filer  (ADCP).  An  empirical  orthogonal  function  (EOF)  decomposition  of  the  fluctuations  was  done,  and  the 
frequency-domain  spectrum  of  the  lowest  empirical  mode  was  presented  in  Fig.  19  of  Ref.  23.  Due  to  an 
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error  in  that  figure,  the  published  power  spectral  density  is  too  large  by  a  factor  of  480  [58].  Corrections  for 
this  error  were  made  before  any  further  processing  of  the  SWARM  spectrum. 

The  spectra  published  for  Barents  [54,  Fig.  9]  and  SWARM  [23,  Fig.  19]  are  both  plotted  as  displacement 
variance  (nr/cpd)  versus  frequency  in  cycles-per-day  (cpd).  These  figures  were  digitized  by  hand  and 
converted  to  units  of  energy  density  per  wavenumber.  First  the  spectra  were  rescaled  by  a  factor  of  86400/27: 
to  convert  the  units  to  m2/  (rad/s). 

Next,  the  spectra  were  converted  from  the  frequency  to  the  wavenumber  domain  under  the  assumption 
that  these  environments  are  dominated  by  the  first  internal-wave  mode,  i.e.,  j  =  1  in  what  follows.  Let  Sj(co) 
with  units  m 2 /  (rad/s)  represent  the  frequency-domain  displacement  spectra.  Then  because  5/ (ft)) dot  = 
Sj(Yh)d-Yh-  the  wavenumber  domain  spectrum  is  Sj(%)  =  Sj((o)vJg(yh),  where  i/,(y/,)  =  d(Oj{%)/dYh  is  the 
modal  group  speed;  Sj(yi,)  now  has  units  m2/  (rad/m).  The  required  group  speeds  were  obtained  by  solving 
the  internal-wave  Sturm-Liouville  problem  [25,  44,  59]  using  the  density  profiles  reported  for  the  experi¬ 
ments’  environments. 

Finally,  the  displacement  spectra  were  converted  to  energy  spectra  by  recognizing  that  Sj{%)  is  the 
spectrum  of  a  variance,  which  can  be  substituted  into  a  wavenumber-domain  transform  of  Eq.  (67), 


Fj(yh)  =  pSj(yh)  fZ° N2(z)dz,  (78) 

Jo 

where  Fj(jh)  is  an  energy  spectrum  with  the  same  units  (J/m2)  /  (rad/m)  as  the  Garrett-Munk  spectrum. 
With  the  above  relationships,  the  frequency-domain  power  spectral  density  presented  in  the  Barents  and 
SWARM  papers  were  digitized  and  converted  to  a  wavenumber-domain  energy  spectrum  for  comparison 
with  the  Garrett-Munk  spectrum. 

5.1.2  TAVEX  Spectra 

The  Transverse  Acoustic  Variability  Experiment  (TAVEX)  [60]  took  place  in  the  East  China  Sea  during 
August  2008.  The  objective  of  the  experiment  was  to  obtain  data  concerning  the  impact  of  dynamic  ocean 
processes  on  the  coherence  of  acoustic  fields  along  a  direction  transverse  to  the  main  path  of  acoustic  prop¬ 
agation.  In  support  of  that,  a  vertical  conductivity,  temperature,  and  depth  (CTD)  chain  with  sensors  at  fixed 
intervals  of  depth  was  towed  along  a  series  of  straight-line  tracks.  Because  the  track  data  were  analyzed  in 
3  km  segments,  obtaining  reliable  estimates  of  low-wavenumber  components  of  the  spectra  are  precluded  in 
this  environment.  However,  these  sensors  were  able  to  provide  measurements  of  high-wavenumber  fluctua¬ 
tions  as  a  function  of  horizontal  position.  From  the  density  fluctuations  the  variance  of  vertical  displacement, 
and  the  spectrum  of  that,  were  obtained  as  follows.  The  relationship  between  density  fluctuations  and  small 
vertical  displacements  due  to  linear  internal  waves  is 

8p(y,z)  =  dp^Z\(y,z),  (79) 

where  8p(y,z)  =  p  (y,  z)  —  p(z)  is  the  density  perturbation,  £(y,z)  is  the  internal- wave  displacement,  and 
dpp(z)/dz  is  the  potential  gradient  of  the  mean  density  profile  p(z)  =  ( p{y,z))y .  Note  that  with  these 
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definitions  both  Sp(y,z)  and  £(y,z)  are  zero  mean.  In  shallow  water  where  pressure  corrections  are  small,  a 
reasonable  approximation  of  buoyancy  frequency  squared  [44,  59]  is 


N2(z)  = 


8  dPp(z) 
p(z )  dz 


(80) 


Equation  (79)  was  solved  for  £( y,z )  and  Fourier  transformed  over  the  transverse  dimension,  y,  to  obtain 
the  power  spectral  density  S(yi,,z)  of  the  displacement  variance  as  a  function  of  horizontal  wavenumber  and 
depth.  As  was  done  for  the  SWARM  data  for  Eq.  (78),  this  was  substituted  into  Eq.  (67),  but  this  time  taking 
depth  variability  into  account: 


F(lh)  =  [  p{z)N2{z)S(yh,z)dz. 
J  o 


(81) 


5.1.3  Spectra  Discussion 

Internal- wave  energy  spectra,  estimated  using  the  above  procedures,  are  shown  in  Fig.  2.  Here  in  the  log- 
log  wavenumber  domain,  the  spectra  lie  along  a  common  power-law  line  with  a  slope  that  is  approximately 
—2.  The  bump  in  the  SWARM  spectrum  between  wavenumbers  1 0  2  and  10  1  rad/m  was  attributed  to 
nonlinear  internal  waves  by  Apel  et  al.  [23].  The  rise  in  the  TAVEX  spectrum  at  the  highest  wavenumbers  is 
likely  aliasing  due  to  contamination  by  fine-scale  features  such  as  spice,  which  were  observed  by  the  CTD 
chain  [61]. 


Fig.  2  —  Empirical  best-fit  Garrett-Munk  spectrum  (black)  with  E 0  =  200  J/m2,  deep-water  Ej  (%), 
and  shallow-water  j *  =  1  with  p  =  4.  For  comparison,  spectra  are  shown  for  TAVEX  (brown), 
SWARM  (magenta),  and  Barents  (blue). 
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A  Garrett-Munk  spectrum  was  empirically  fit  to  the  experimental  data.  The  parameters  that  best  describe 
the  spectrum  include  a  large  energy  density  of  £0  =  200  J/m2,  which  is  seemingly  more  appropriate  for 
shelf-break  regions  with  water  depths  on  the  order  of  200  m,  than  shallow  shelf  areas  where  the  energy  is 
expected  to  be  an  order  of  magnitude  less.  The  spectral  slope  is  also  odd,  being  best  fit  with  the  deep-water 
spectrum  for  Ej{jh )  given  by  Eq.  (73).  The  choice  of  shallow-water  coefficients,  j,  =  1  and  p  =  4,  merely 
emphasizes  the  first  mode  above  all  others. 

It  must  be  noted  here  that  both  the  SWARM  and  TAVEX  experiments,  though  situated  upslope,  were 
in  the  vicinity  of  shelf  breaks  at  200  m  depth  where  internal  waves  are  known  to  be  generated  [52,  62,  63], 
and  the  water  depth  across  the  Barents  site  spanned  100-400  m.  So  the  parameters,  although  unusual,  are 
not  entirely  unjustifiable.  And  while  the  parameters  of  this  fit  may  be  debatable  on  physical  grounds,  the 
resulting  spectrum  does  mimic  the  data  closely  enough  to  provide  a  reasonable  simulation  of  the  acoustic¬ 
mode  fluctuation  spectrum,  Eq.  (75),  that  is  required  by  the  phase-structure  function,  Eqs.  (76)  and  (77). 

5.2  Coherence  Discussion 

Prior  to  discussion,  a  few  words  need  to  be  said  about  acoustic  and  internal-wave  mode  functions.  To 
determine  the  acoustic  coherence  function  it  is  necessary  to  evaluate  the  modal  interaction  described  by 
Eq.  (57),  which  requires  suitable  acoustic  and  internal-wave  modes.  Because  the  TAVEX  environment  was 
recently  measured  by  the  author,  it  was  chosen  to  provide  those  modes.  Acoustic  modes  were  computed  by 
the  Kraken  program  of  Porter  [64],  and  internal- wave  modes  were  computed  using  a  program  by  Bell,  Jr. 
[65], 

Computation  of  the  acoustic  coherence  function  also  requires  the  partition  of  acoustic  energy  among 
modes  as  prescribed  by  Eq.  (21).  The  necessary  mode  amplitudes  were  obtained  by  projecting  a  complex 
acoustic  field  from  the  RAM  parabolic  equation  model  of  Collins  [66-68]  onto  the  acoustic  mode  functions 
from  Kraken.  Figure  3  shows  the  partition  of  intensities  among  modes  at  500  Hz.  Clearly  the  lowest  acoustic 
modes  have  the  most  weight,  so  their  properties  tend  to  dominate  the  combined  acoustic  coherence. 

In  comments  following  Eq.  (66),  it  was  pointed  out  that  the  phase-structure  function  increases  with  mode 
number.  This  behavior  is  seen  in  Fig.  4  where  the  amplitude  of  the  phase-structure  function  increases  much 
more  rapidly  for  mode  7  (black)  than  for  mode  1  (blue),  which  remains  flat-lined.  Also  apparent  in  this 
figure  is  the  approach  to  saturation  at  a  range  of  5  km  with  decaying  undulations  caused  by  the  sine  function 
in  Eq.  (46).  Obviously,  the  modal  coherence  functions  must  mirror  this  behavior  (Fig.  5)  because  they  are 
the  negative  exponentials  of  the  phase-structure  functions,  cf.  Eq.  (20). 

The  total  acoustic  coherence  shown  in  Fig.  6  is  the  weighted  sum,  Eq.  (21),  of  the  modal  coherences. 
In  this  figure  it  is  seen  that  the  coherence  falls  off  more  rapidly  with  transverse  separation  at  500  Hz  than 
at  300  Hz.  This  observation  is  in  line  with  the  comments  following  Eq.  (66),  where  such  behavior  was 
predicted.  In  comparing  the  500  Hz  total  coherence  (blue)  to  the  modal  coherences  in  Fig.  5  it  is  seen  that 
the  total  coherence,  as  the  weighted  sum  of  all  modes,  behaves  in  this  case  much  like  mode  4  (cyan). 

The  power-law  behavior  of  the  coherence,  n,  discussed  at  the  end  of  Sec.  2.1,  can  be  obtained  by  twice 
taking  the  logarithm  of  Eq.  (9)  and  substituting  Eq.  (10), 

Ap(y)  =  ln(— 21n(Rp(Ay)))  =lnDP(Ay), 

=  nlnAy  —  n\nL. 


(82) 

(83) 
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Fig.  3  —  Partition  of  acoustic  intensity  among  modes  for  500  Hz  signals 


Fig.  4  —  Modal  phase-structure  functions  for  500  Hz  acoustic  modes  1-7 
colored  in  order  blue  to  black 
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Fig.  5  —  Modal  coherence  functions  for  500  Hz  acoustic  modes  1-7  colored  in  order  blue  to  black 


Fig.  6  —  Phone  coherence  at  300  Hz  (black)  and  500  Hz  (blue) 
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The  resulting  theoretical  power-law  behavior,  shown  in  Fig.  7,  has  a  slope  of  n  =  —2  for  small  separations, 
as  expected  from  the  properties  of  the  sine  function  in  Eq.  (46).  But  this  slope  flattens  out  as  the  coherence 
saturates  to  the  scattering  intensity  at  large  separations.  The  simple  power  law  described  by  Eq.  (10)  only 
holds  for  small  separations. 


Fig.  7  —  Log-log  coherence  power  law  A/>(ln(y))  from  Eq.  (83)  shown  for 
300  Hz  (black)  and  500  Hz  (blue) 


The  coherence  presented  in  these  numerical  evaluations  of  path-integral  theory  decays  much  more 
slowly  than  is  found  in  experimental  measurements.  These  results  are  clearly  at  odds  with  a  significant 
body  of  the  literature.  Carey  [29]  reported  a  comprehensive  list  of  coherence  lengths  in  various  seas  in  the 
135-600  Hz  band,  with  deep-water  lengths  on  the  order  of  1 00 A  over  500  km,  and  shallow-water  coherence 
about  30 A  over  45  km.  In  the  Levantine  Sea,  Carey  [69]  found  coherence  lengths  on  the  order  of  543,  231, 
and  116  m  at  frequencies  of  58,  175,  and  348  Hz,  respectively.  For  observations  of  transverse  coherence 
at  a  range  of  40  km  in  the  Strait  of  Korea,  Carey  et  al.  [30]  found  coherence  lengths  of  1 12  m  at  400  Hz. 
For  deep-water  data  obtained  in  the  Pacific,  Andrew  et  al.  [15]  observed  coherence  lengths  on  the  order 
of  410-528  m  at  ranges  between  2000  and  3000  km.  These  results  were  obtained  assuming  a  power-law 
phase-structure  function  like  Eq.  (10)  with  n  =  3/2  that  was  based  on  the  work  of  Flatte  and  Stoughton  [26]. 


6.  SUMMARY 

The  mutual-coherence  function  was  found  in  Eq.  (8)  to  be  described  by  a  normalized  amplitude  cor¬ 
relation  function  multiplied  by  the  phase-structure  function.  In  shallow  water  where  the  acoustic  field  is 
naturally  described  by  a  linear  combination  of  normal  modes,  the  mutual-coherence  function  was  shown  to 
be  a  sum  of  modal  coherence  functions,  each  of  which  is  governed  by  a  modal  phase-structure  function. 
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Formulation  of  the  transverse  coherence  problem  using  mode  amplitudes  propagating  in  the  horizontal 
plane  is  advantageous  because  the  ray  paths  are  nearly  straight  lines  with  no  interaction  with  lateral  bound¬ 
aries.  These  conditions  justify  using  the  parabolic  and  Markov  approximations  required  by  the  path-integral 
approach.  Conversely,  the  modal  approach  breaks  down  near  caustics  that  form  in  the  vicinity  where  modes 
become  evanescent. 

For  isotropic  internal-wave  fields,  fluctuations  with  wavevectors  parallel  to  the  acoustic  path  were  found 
in  Sec.  3.3  to  produce  a  delta  function,  5(y/,cos(0)),  that  selects  only  the  9  =  ±n/2  components  of  the 
internal- wave  spectrum,  thus  showing  that  the  transverse  acoustic  coherence  depends  solely  on  the  transverse 
part  of  the  internal-wave  spectrum. 

The  impact  of  internal  waves  on  acoustic  modes  is  proportional  to  the  internal-wave  energy  £0,  and 
was  found  in  Eq.  (57)  to  depend  on  how  the  shape  of  the  acoustic  mode  function  aligns  with  the  internal- 
wave  excitation.  This  interaction  is  strongest  for  the  lowest  internal-wave  modes,  which  have  maxima  in 
regions  of  high  buoyancy,  and  is  weakest  for  high-order  internal-wave  modes  that  oscillate  rapidly  at  the 
same  depth.  Likewise,  acoustic  modes  with  maxima  in  the  vicinity  of  high  buoyancy  suffer  more  severe 
coherence  loss  than  modes  with  a  null  near  the  gradient.  The  highest  acoustic  modes  are  also  subject  to 
rapid  loss  of  coherence  because  of  the  inverse  dependence  of  the  phase-structure  function  on  the  horizontal 
acoustic  wavenumbers,  as  shown  in  Eq.  (66). 

Although  the  theory  presented  herein  is  perhaps  less  than  fully  satisfying  because  of  the  simplifications 
required  to  make  progress  (adiabatic  modes,  homogeneous  isotropic  internal-wave  spectra,  the  Markov  ap¬ 
proximation,  and  the  Garrett-Munk  spectrum),  it  is  nevertheless  a  complete  theory  that  quantitatively  relates 
the  transverse  acoustic  coherence  function  to  the  underlying  internal-wave  field.  As  such,  it  provides  predic¬ 
tive  capabilities  that  illuminate  several  aspects  of  the  relationship.  Chief  among  these  is  the  result  expressed 
by  Eq.  (46)  that  the  transverse  acoustic  coherence  is  determined  by  transverse  fluctuations  in  the  underlying 
sound-speed  field,  as  expressed  by  an  integral  over  the  fluctuation  spectrum  weighted  by  a  sine  function. 
It  is  the  sine  function  that  provides  the  correct  transition  between  the  phase-structure  function  at  small  hy¬ 
drophone  separations  and  the  scattering  strength  at  large  separations. 

Results  show  the  coherence  predicted  by  path-integral  theory  as  applied  to  the  transverse  scattering  of 
normal  modes  appears  to  decay  too  slowly  with  hydrophone  separation  to  describe  the  observed  ocean. 
In  contradistinction  to  experimental  measurements,  saturation  effects  described  by  the  theory  only  become 
significant  at  very  large  separations  spanning  several  kilometers. 
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Appendix  A 


EIGENVALUE  PERTURBATION  EXPANSION 


Following  Rajan  et  al.  [Al],  a  perturbation  technique  introduced  by  Schrodinger  [A2]  and  discussed 
in  several  standard  references  [A3,  A 4]  is  used  to  relate  stochastic  fluctuations  in  the  eigenvalues  of  the 
Helmholtz  equation  to  those  of  the  underlying  index  of  refraction.  First,  perturbation  expansions  for  the 
eigenvalues  Km  are  introduced  using  Eq.  (25),  which  has  deterministic  part  km  (1  +  um{x.y))  and  small 
stochastic  part  km}im(x.y).  Likewise,  the  total  wavenumber  k  is  expanded  by  letting 


k(x,y,z)  =k(  1  +v(x,y;z)  +  v(x,y,z))  ■ 


(Al) 


Here  1  +v(z)  ~  A  (  I )  is  the  deterministic  wavenumber  profile  and  v(z)  is  a  small  stochastic  fluctuation  such 
that  viz)  ~  im)  AC  1.  The  eigenfunctions  are  represented  as 


T 'm(x,y,z)  =  xFm(z)  +  HmYm(x,y',z) 


(A2) 


where  x¥m(z)  are  the  deterministic  eigenfunctions,  and  the  /i„,  y/)(, (.r,y;z)  are  small  stochastic  functions. 
A  solution  for  the  eigenvalue  perturbations  is  obtained  by  rewriting  Eq.  (14)  in  operator  form  as 


[H (z)  -  Wm]  vU/„  (z)  =  0,  where 


(A3) 


and 


(A4) 


W  —  k2 

yvm  —  /v/7r 


(A5) 


Using  the  above  perturbation  expansions  and  dropping  quadratic  products,  the  problem  can  be  split  into  a 
large  deterministic  part  (0)  and  a  small  stochastic  part  (1): 


H{z) 

Wm 

%n(z) 


(A6) 

(A7) 

(AS) 
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where  the  various  parts  are  defined  as: 


//(0)(z) 

-  piz)Sz{p(z)Sz)+kHl  +  Mz))’ 

(A9) 

=  2k2v(z), 

(A10) 

w<?> 

=  (  1  “1“  2wm), 

(All) 

=  2kliim, 

(A  12) 

^n\z) 

=  %n(z), 

(A  13) 

^(z) 

=  tym  (^)  • 

(A  14) 

Following  established  procedures  [A3,  A4],  first-order  terms  are  collected  and  found  to  satisfy 


'Pm'1  (z)/p(z)dz  =  o. 


(A  15) 
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Appendix  B 

LIST  OF  SYMBOLS 


<>t 

ensemble  average  over  time 

<>y 

ensemble  average  over  space 

<>e 

ensemble  average  over  environments 

X 

range  from  source 

e 

transverse  polar  coordinate 

y  =  re 

transverse  Cartesian  coordinate 

z 

depth  below  ocean  surface 

x0 

receiver  range 

Zo 

receiver  depth 

Ay 

transverse  hydrophone  separation 

co  =  2k  f 

acoustic  frequency 

U  =  (r,0,z) 

generalized  coordinate  vector 

c(  u) 

sound  speed 

c 

mean  sound-speed  scale 

K  u) 

acoustic  wavenumber 

P(  u) 

complex  acoustic  pressure 

P(y) 

complex  acoustic  phase 

aj 

acoustic  pressure  intensity 

< 

acoustic  mode  intensity 

L 

empirical  coherence  length 

n 

empirical  coherence  power-law 

scattering  strength 

8(j) 

phase  fluctuation 

X 

arbitrary  Cartesian  coordinates 

s 

arc  length  along  unperturbed  path 

k(x) 

unperturbed  wavenumber 

D(x  i,x2) 

generalized  phase-structure  function 

DP  (Ay) 

transverse  acoustic  phase-structure  function 

Dm(Ay ) 

transverse  modal  phase-structure  function 

Ap(Ay) 

logarithmic  phase  structure  In  Dp(Ay) 

RP(Ay) 

acoustic  pressure  coherence 

RAm(Ay) 

acoustic  mode  coherence 

RP(Ay) 

acoustic  phase  coherence 

RA,„(Ay ) 

modal  phase  coherence 

Hm 

acoustic  mode  structure  weight 
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Qo 

x¥m(r,d-,z) 
Am  (p  0  ) 
Am  (v) 

Am 

km(r,d) 

Km(x,y) 


Um(x,y) 
Hm(x,y) 
v(x,y;z ) 
v(x,y,z) 

Pm 

P 

8p 

¥m(x,y) 

le 

Dp 

Um(x,yp(x)) 
a,  j8 

A, 

ya(x) 

Aap(r) 

Y 

Ya,Yp 

Ay 

Yh 

Yx,Yy 

Am  (  Y*  i  %  ) 
da/dz 

dp/dz 

fc 

N(x,y,z ) 

N 

-  Ain.j  ( Yh ) 

C(*,:y;z,0 
Wj(x,y;z,Yh ) 


acoustic  source  strength 

acoustic  mode  eigenfunction 

acoustic  mode  amplitude 

mode  amplitude  phase 

reduced  mode  amplitude 

acoustic  horizontal  wavenumber  (eigenvalue) 

reduced  acoustic  horizontal  wavenumber 

mean  reduced  horizontal  wavenumber 

horizontal  wavelength  of  acoustic  mode 

deterministic  eigenvalue  perturbations 

stochastic  eigenvalue  perturbations 

deterministic  sound-speed  perturbations 

stochastic  sound-speed  perturbations 

two-point  correlation  of 

mean  density  scale 

zero-mean  density  fluctuation 

parabolic  approximation  of  mode  amplitude 

environmental  fluctuation  scale  length 

integration  over  all  paths 

deterministic  paid  of  path  integral 

principal  ray  labels 

mean  range 

range  separation 

transverse  position  of  rays 

transverse  separation  of  rays 

transverse  position  along  array 

transverse  position  of  rays  at  array 

transverse  separation  at  array 

horizontal  wavenumber  of  environmental  fluctuations 

Cartesian  wavenumbers  of  environmental  fluctuations 

environmental  fluctuation  spectrum 

adiabatic  gradient 

potential  gradient 

local  Coriolis  frequency 

internal-wave  frequency 

buoyancy  frequency 

depth- averaged  buoyancy  frequency 

buoyancy-modal  interaction  integral 

internal-wave  displacement 

internal-wave  eigenfunctions 

internal-wave  spectral  amplitude 
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Gj(yh) 

g 

Ej(Yh) 

Eo 

Ej  ( Yh ) 

Qj 

Mj 

M 

j* 

P 

Yj 

Sj(co) 
Ej(Yh ) 


internal-wave  power  spectra 
internal-wave  structure-function  integral 
internal-wave  energy  per  unit  area 
Garrett-Munk  spectrum 
Garrett-Munk  energy  scale 
normalized  Garrett-Munk  spectrum 
internal-wave  spectral  normalization  constant 
internal-wave  mode  energy  partition 
internal-wave  mode  normalization  constant 
Garrett-Munk  reference  mode 
Garrett-Munk  mode  exponent 
Garrett-Munk  spectral  peak 

experimental  internal- wave  spectrum  (m2/  (rad/s)) 
experimental  internal- wave  spectrum  ((J/m2)  /  (rad/m)) 


